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Abstract. - We investigate the transport of energy, magnetization, etc. in several finite 
one- dimensional (ID) quantum systems only by solving the corresponding time-dependent 
Schrodinger equation. We explicitly renounce on any other transport-analysis technique. Vary- 
ing model parameters we find a sharp transition from non-normal to normal transport and a 
transition from integrability to chaos, i.e., from Poissonian to Wigner-like level statistics. These 
transitions always appear in conjunction with each other. We investigate some rather abstract 
"design models" and a (locally perturbed) Heisenberg spin chain. 



The transport behavior of one-dimensional (ID) systems has intensively been investigated 
for several decades, as well in the context of classical mechanics as in the context of quan- 
tum mechanics [1-16]. Nevertheless, the precise conditions under which normal transport 
occurs, i.e., under which there is neither ballistic transport nor localization but normal spatial 
diffusion, are still not known [17]. 

In the classical domain it seems to be largely accepted that normal transport (in any dimen- 
sion) requires the chaotic dynamics of a non-integrable system whereas non-normal transport 
is typical for the regular dynamics of (completely) integrable systems, see [3]. However, there 
have also been successful attempts to observe normal transport in the absence of exponential 
instability, the latter being a basic feature of (deterministic) chaos [4]. In the quantum do- 
main there are only very few examples which can be reliably shown to exhibit normal, diffusive 
transport at all [5-7]. But, also in this field, it has been argued that non-normal transport 
is related to the macroscopic number of conserved quantities which characterize integrable 
systems [10-13]. Moreover, recent numerical computations for spin chains have led to the 
assumption that normal transport might strictly depend on quantum chaos [1,2]. 

Although this assumption is plausible, it has not been proved yet [15]. Moreover, almost all 
computations are either based on special models of reservoirs which might effect the transport 
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behavior, or they rely on the Kubo formula. The latter has originally been derived for field- 
driven electrical conductance and its validity for diffusive transport phenomena such as, e.g., 
thermal conductance is still under dispute, see [7-9,18]. We refer to [19] where the rather lim- 
ited validity of the Kubo formula for thermal conductance has explicitly been demonstrated. 
Especially considerations which are restricted to the analysis of the so-called "Drude-weight" 
are not sufficient to determine the transport behavior. Thus the main intent of the letter 
at hand is to examine the relation between transport behavior and quantum chaos without 
modelling, e.g., external heat baths at different temperatures or using the Kubo formula. 

This letter is structured as follows: First of all we briefly comment on the theory of quan- 
tum chaos, mainly on the nearest neighbor level spacing distribution (NNLSD). Thereafter 
we introduce two measures for deviations: i) a measure for the deviation of a given system 
from the fully chaotic Wigner NNLSD (xw), h) a measure for the deviation of a given sys- 
tem from fully diffusive behavior (D). Different finite ID models are considered and for each 
model a parameter that drives an "integrable to chaotic transition" is varied. The corre- 
sponding Schrodinger equations are solved. This allows to plot both measures (xw-> D) over 
the respective parameter in order to reveal correlations. This reveals that, at least in our 
models, the "integrable to chaotic" induces a "non-normal to normal transport" transition. 
We investigate some "design models" featuring random interactions and a (locally perturbed) 
Heisenberg spin chain. We close with a summary and a conclusion. 

The theory of quantum chaos is principally concerned with the level statistics of quantum 
systems which possess a classical limit [20-26]. A commonly used statistical measure is the 
nearest neighbor level spacing distribution (NNLSD), P(s), where P(s)ds is the probability 
that the distance, s, between two adjacent eigenvalues lies in the interval [s, s + ds]. Typically, 
P(s) is well described by a Wigner distribution Pw(s) = ns/2 exp(— 7rs 2 /4), when the clas- 
sical limit is chaotic, and by a Poissonian distribution Pp(s) = exp(— s), when the classical 
limit is regular, i.e., (completely) integrable [20,21]. Interestingly, P(s) can differ from these 
distributions [26]. Two things are crucial for the computation of the NNLSD: First of all one 
has to select a subspace consisting of states from a single symmetry class. Thereafter one 
has to unfold the subspace's spectrum, such that the local average of s equals one (s = 1) 
everywhere in the spectrum. A detailed description of the unfolding procedure can be found 
in [20,21]. In order to compare the resulting NNLSD, given as a normalized histogram with 
L bins, with the above distributions we define the measure 

and xp: respectively, where P M is the probability that s lies inside the /x'th bin of the histogram. 
Pw,(i an d Pp,n are t ne probabilities according to Pw(s) and Pp(s), respectively. 

In this letter we consider chain-like quantum systems which may be described by Hamil- 
tonians of the form 

N N-l 

& = E k + E *w+i . (2) 
f-i—i i 

where denotes the local Hamiltonian of some subunit /i and v^.^+i the interaction between 
neighboring subunits, N is the total number of subunits. The model is primarily intended 
to investigate energy transport in a chain of coupled quantum systems (like, e.g., molecules) 
but may also be viewed as a Hubbard- type model for particles on a lattice (see below). 
Furthermore, most (almost) periodic systems should allow for a description according to J2J, 
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Fig. 1 - A chain of N identical subunits: Each subunit features a nondegenerate ground state, a 
wide energy gap (AE) and a comparatively narrow energy band (Se) which contains n, energetic 
equidistant, states. 



as will be illustrated below with the example of a spin chain and is explained in detail in [19]. 
Since we intend to identify the diffusion of, e.g., energy, we have to introduce a measure for 
the energy density or the local energy. We define the local energy operator at site fi simply as 
hp. Of course, this definition neglects the energy contained in the interaction and eventually 
implies a weak coupling limit. 

If the transport behavior of, e.g., energy was perfectly diffusive, then the local energies 
Efj,(t) — (^(t^h^ij}^)) (\ip(t)) being the full system's wavefunction) would obey the following 
set of equations for, e.g., a chain- like system 

Ex =r ] (E 2 -E 1 ) , 

E^ = v [(E^ +1 - E„) - (E„ - , (3) 

En = V (En -i - E N ) ■ 

This may be viewed as a discrete form of the diffusion equation pe = r/ApE, Pe being the 
energy density. Diffusive behavior of any other quantity may be defined by a respective set of 
equations. 

A suitable measure for the deviation of the (numerically) exact time evolution of the 
local energies as obtained from the Schrodinger equation from the fully diffusive dynamics as 
generated by J3J is given by 

D2 = N^J2J Q [E™ ma \t)-E^(t)] 2 dt (4) 

with t = 1/t). In the letter at hand we determine rj through fitting the dynamics generated 
by © to the (numerically) exact time evolution of the local energies. For a possible analytic 
computation of 77, see, e.g., [6]. 

Let us now investigate two abstract examples ("design models") of The first model 
consists of N identical subunits: Each subunit features a nondegenerate ground state, a wide 
energy gap (AE) and a comparatively narrow energy band (Se) which contains n, energetic 
equidistant, states, cf. Fig. ^ The next-neighbor interaction is defined as 

n 

where h.c. is the hermitian conjugate of the previous sum. corresponds to an upwards 
transition of the /i'th subunit from its ground state to the i'th state of its band and ff^ i 
corresponds to a downwards transition, respectively, are randomly distributed complex 
numbers which are normalized to - =1 \vij\ 2 /n 2 = 1, such that A sets the total interaction 
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Fig. 2 - (a) Average deviation of the (numerically) exact time evolutions of the local energies from 
perfectly diffusive behavior, (b) Deviation of the level statistics from the Wigner distribution for the 
model pictured in Fig. Q The model parameters N = 2, AE — 10 and Se = 0.5 are fixed, but n 
runs from 125 (+) over 250 (x) and 500 (*) till 1000 (□). Evidently, Wigner-like level statistics and 
normal transport are correlated. 



strength. Since does not change with fj,, there is no disorder. Obviously, we can restrict 
our analysis to the one-excitation subspace (the space where one subsystem is excited and 
all others are in their ground state), if the initial state ip(0) belongs to this subspace. The 
subspace's dimension is N ■ n and grows linearly rather than exponentially with N. 

The model may be also illustrated as a "single-particle multi-channel quantum wire" , as 
already mentioned above, with random hoppings but without disorder. Therefore, unlike the 
single-particle Hubbard models with randomness in, e.g., [27], no localization occurs. Thus 
all deviations from normal transport which are discussed below (see Fig. are deviations 
towards a ballistic- type behavior, not towards localization. 

For simplicity, we start with N — 2, i.e., only two subunits. Figure^ shows D 2 versus 
A, averaged for 200 adequately restricted random pure initial states i/j(0), cf. The model 
parameters AE — 10 and Se = 0.5 are fixed, but n runs from 125 to 1000. FigureEJ} shows Xw 
versus A, that is, the deviation of the corresponding NNLSD from the Wigner distribution, cf. 
|QJ. Expectedly, for too small A (< 10 -5 ) the subunits are almost uncoupled, such that the 
time evolution of the local energies differs from J3J and the NNLSD deviates from the Wigner 
distribution. Evidently, © does not apply for too large A (> 10~ 2 ), although the NNLSD is 
clearly Wigner-like. Obviously, a Wigner-like NNLSD is not a sufficient condition for normal 
transport. Nevertheless, for all n one observation is striking: The minimum of D 2 lies exactly 
at the position where the minimum of Xw 1S reached. This observation which turns out to 
apply to all our models is our main result. Thus, a Wigner-like NNLSD might be a necessary 
condition for normal transport. 

We have additionally checked chains up to N = 15 and n = 500. Furthermore, we have 
varied details of the model, e.g., the band's level distribution. We shortly summarize that the 
results do not significantly differ from the results of the case N = 2. The interested reader is 
referred to [6, 19] where, e.g., the finite size scaling is discussed in detail. 

However, since real physical systems do not typically feature the gapped local spectra of the 
above model, we consider another example of @- The model is pictured in Fig.|21and consists 
of N identical subunits: Each subunit features n eigenstates which are randomly, but uniformly 
distributed within an energy interval AE. As before u^.^+i is a randomly chosen complex 
matrix, but now without restriction to any subspace (no "particle- number conservation"). 
Nevertheless, is supposed to be given by v^^+i = At) with Tr(v 2 )/n 2 = 1. Again, for 

each A we average over an adequate set of 200 random pure initial states. For the case N = 2 
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Fig. 3 - A chain of AT identical subunits: Each subunit features n states which are randomly, but 
uniformly distributed within an energy interval AE. 

the results do not qualitatively differ from the results in Fig. and[2b, e.g., for the model 
parameters n = 60 and AE = 10 the minimum of D 2 and Xw ^ es at A = 0.0005. 

What about larger TV? Since for this system class the dimension of the relevant Hilbert 
space is n , we are forced to decrease n. But for smaller n deviations from the fully diffusive 
behavior typically increase, cf. Fig.|2Ji. Nevertheless, in principle we find the above result also 
confirmed for N = 3 and n = 20. Remarkably, even for N — 6 and n — 4 the time evolution of 
the local energies is in tolerably good agreement with Q , again for a parameter regime where 
the NNLSD is Wigner-like. (We do not display all those data here, since they essentially look 
like Fig. El) 

However, for a chain of two level subunits, that is, for a chain of "spins" relaxation and 
local fluctuations are indistinguishable on the scale of a single subunit, i.e., on this scale 
J2J does not apply. This changes if the scale is changed, i.e., if various neighboring spins 
(including their mutual interactions) are grouped together to form a subunit as addressed by 
However, since a spin chain with random next neighbor interaction may always be viewed 
as a "mixture" of different wellknown spin models, it is more meaningful to consider directly a 
single spin model, e.g., the Heisenberg model. Furthermore, other than our "design models", 
those spin models possess a classical counterpart which is unambiguously either integrable 
or not [21,28]. Thus, in the following we consider a Heisenberg spin chain in an external 
magnetic field B. Concretely, the Hamiltonian reads 

JV N-l 

^ = 2 E( S + + 4 E ^ ' fr/x+i . N even, (6) 

where are the standard Pauli operators and are local variations from B. Furthermore, 
-B M are chosen as Gaussian distributed random numbers with ( B^ ) = and {B^B U ) = e 2 . 
The parameter A sets the coupling strength, but the Heisenberg interaction is not normalized 
to 1, unlike JSJ. As mentioned above we operationally divide the chain into only two subunits, 
namely the first and second half. Note that due to the local fields B^ the spectra of those 
halves may not be identical. 

Obviously, © is invariant under rotations around the z-axis. In order to compute the 
NNLSD we choose the subspace with M = 0, where M is the quantum number with respect 
to S z = Ylu=i <^u/2j the generator of these rotations. For a model as given by © with N = 12 
the dimension of this subspace is d — 924. Figure EJi shows Xw versus e i averaged for 100 
sequences B^, that is, the average deviation of the NNLSD from the Wigner distribution, 
and Xpi respectively, cf. QJ. Since the system is integrable for e = 0, a Poisson-like NNLSD 
is obtained. When e increases from zero, the system undergoes a transition to chaos and 
consequently the NNLSD becomes Wigner-like. The minimum of Xw (the maximum of Xp) 
is reached at e « A/4. When e further increases and becomes larger than A, the system 
becomes localized and accordingly a Poisson-like NNLSD reappears, see [22-24], too. 
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Fig. 4 - (a) Average deviation of the level statistics from the Wigner distribution (+) and the Poisson 
distribution (x) for the spin model corresponding to (JSJ. (b) Time evolution of the first half's 
local magnetization for e/A = (dashed curve), 0.25 (dotted curve) and 1 (solid curve). Evidently, 
exponential decay indicating normal transport occurs only for e/A = 0.25, the parameter for which 
the level statistics are most Wigner-like. 



Since neither the level statistics nor the dynamics within the mentioned subspace depend 
on the constant field B, we set B = 0. Here we analyze the transport of magnetization which 
is, just like energy, a conserved quantity in this model. According to our above partition 
scheme we define two local magnetizations: 

Mi:=l(m\j2^M^ and M *-=\m)\Y,°i\m)- (?) 

£1=1 fJ,=7 

If the transport behavior of magnetization was diffusive on this scale, the M's should exhibit a 
dynamics as generated by a direct analogue to i.e., they should simply relax exponentially 
to equilibrium. In Fig. ^3 Mi is displayed for the cases e/A = (dashed), 0.25 (dotted) 
and 1 (solid). The chosen initial state i/j(0) is the only state with Mi = 3 and Mi = —3. 
Obviously, for the case e/A = the transport is non- normal, the "bouncing" behavior of the 
magnetization could rather be interpreted as a hint for ballistic transport. This issue has 
been discussed very controversially in the literature [12-16]. So far, we also have no definite 
conclusion. For e/A = 1 almost no transport is observable, the magnetization seems to be 
stuck, i.e., localized. But at the minimum of Xw (^ ne maximum of xV), e/A = 0.25, there 
is a tolerably good agreement with JSJ: the first half's local magnetization decays almost 
exponentially from the initial value Mi = 3 to the equilibrium value Mi = 0, as expected for 
normal transport. 

What is to be expected for longer chains? Regular transport and localization scenarios 
should be unaffected in the thermodynamic limit as the investigations of the first system class 
and theories in [6,19,29] suggest. The ballistic-like transport may possibly become normal on 
a larger scale. 

Of course, our Heisenberg model can be mapped on a ID twelve-site Hubbard model 
of interacting spinless fermions with disorder. Such systems have recently been studdicd 
in [30] with the result that the interaction may lead to finite conductivity (above some critical 
temperature) for systems that would be localized (in real space) otherwise. Although our 
states are far from being thermal, our results are in accord with those findings: If we remove 
the interactions (the <7*<7* +1 -terrn) for the normal transport case (e/A = 0.25) we find that 
magnetization essentially stays in the first half, i.e., localization occurs. Consistently, the 
NNLSD becomes Poissonian again. 
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Let us finally summarize and conclude: We investigated several finite one-dimensional 
quantum systems. We introduced definitions for perfectly diffusive, normal transport behavior 
and for the level statistics expected from a perfectly chaotic system (Wigner statistics). By 
numerically solving the corresponding Schrodinger equations we got the exact energy spectra 
as well as the exact dynamics of our systems. With those data we computed the deviation of 
our model's level statistics from a perfectly chaotic system's level statistics and the deviation 
of our model's exact transport dynamics from a perfectly diffusive transport dynamics. From 
comparing those deviations we found that models featuring nearly perfect diffusive transport 
always feature nearly perfect chaotic level statistics, while the inverse is not true. This result 
is eye-catching and eventually leads to our conclusion that a Wigner-like NNLSD might be a 
necessary, but not sufficient condition for normal transport. 

* * * 

We thank H.-J. Schmidt for fruitful discussions and gratefully acknowledge financial sup- 
port by the Deutsche Forschungsgesellschaft. 
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